Lie-Nambu and beyond 

Marek Czachor 

Katedra Fizyki Teoretycznej i Metod Matematycznych 

Pohtechnika Gdanska, ul. Narutowicza 11/12, 80-952 Gdansk, Poland 

E-mail: mczachor@sunrise.pg.gda.pl 
*** 

An extended version of the talk given in July 1997 in Peyresq, France 



OS 

t— i 

> 
o 

(N 



> 

o 



a: 
i 

^— > 

G ■ 

3 : 



Linear quantum mechanics can be regarded as a particular 
example of a nonlinear Nambu-type theory. Some elements 
of this approach are presented. 



I. INTRODUCTION 

At the moment there is no single experimental data 
suggesting that states of a quantum system can evolve 
in a fundamentally nonlinear way. On the other hand all 
impossibility theorems stating that such a nonlinearity is 
in principle impossible have not survived a detailed anal- 
ysis. It is therefore possible that the status of quantum 
linearity is similar to this of geometrical linearity from 
before the invention of general relativity. 

The multiple-bracket dynamics described in this paper 
arose from a search for a consistent embedding of linear 
quantum mechanics into a more general theory where the 
assumption of linearity could be dropped. The formalism 
is essentially based on density matrices and not on wave 
functions. A density matrix plays here a role of a funda- 
mental field and should not be regarded as a mixture of 
classical and quantum probabilities. A departure point 
for the discussed generalization is the observation that 
density matrices of ordinary quantum mechanics satisfy 
an equation of a Lie-Nambu type. 

The layout of the paper is as follows. Sec. II re- 
lates the work to the earlier ones by Nambu (1973), 
and Bialynicki-Birula and Morrison (1991). Some for- 
mal tools are introduced in Sec. Ill, IV and V. In Sec. 
VI a (2n + l)-bracket is introduced and some of its gen- 
eral properties are proved. The bracket differs from the 
so-called generalized Nambu, generalized Poisson, or gen- 
eralized Nambu-Poisson brackets discussed in the liter- 
ature (Bayen and Flato, 1975; Cohen, 1975; Takhta- 
jan, 1994; Gautheron, 1996; Chatterjee and Takhtajan, 
1996; de Azcarraga et al., 1996; Ditto and Flato, 1997; 
Ibanez et al., 1997; Ditto et al., 1997; Kanatchikov, 
1997). The particular case of a 3-bracket dynamics is 
discussed in Sec. VII where the notion of a Lie-Nambu 
duality is introduced and properties of solutions of the 3- 
bracket equations are discussed. The 5-bracket dynamics 
is briefly discussed in Sec. VIII. Sec. IX is devoted to 
the question of iV-particle extensions of a nonlinear Lie- 
Poisson dynamics. The notion of complete separability 
is discussed and examples of completely separable equa- 
tions are given. The results of this section contradict the 



popular belief that all nonlinear extensions of quantum 
mechanics lead to faster-than-light phenomena. Also the 
question of complete positivity of solutions is shortly dis- 
cussed in this Section. Sec. X is devoted to the problem 
of separability of the dual Poisson dynamics and it is 
shown that a surprising nonlocal phenomenon occurs. In 
Sec. XI a possible link between our formalism and the 
problem of quantization of classical Nambu dynamics is 
discussed. 



II. LIE-POISSON AS A LIE-NAMBU 

The origin of this work goes back to two papers where, 
in completely different contexts, a notion of a triple 
bracket was introduced. 



A. Nambu 1973: Euler equations 



The Euler equations for a rotating rigid body are 



Jk — £kbc{<h/ h)Jc 

dJk3H_ _ 
' dJ a dJb 
dJ k dH dS 



Ja^ b 



— ^abc 



dJ a dJb dJ c 



{Jk,H} 



(1) 
(2) 

(3) 



Here J is an angular momentum, 1^ a component of a 

J 2 J 2 t/ 2 

moment of inertia, H — ^jj- + ^ + jj- rotational energy, 
and S = \J 2 ■ The totally antisymmetric tensor e a b c can 
be regarded either as a 3-dimensional volume form, or 
as structure constants of so(3). The Lie algebra so(3) 
enters the equations also via S since J 2 is a second order 
Casimir invariant of this algebra. The form (||) defines 
a Poisson bracket. The triple bracket defined by (||) is 
nowadays called the Nambu bracket and was introduced 
in Nambu (1973). 

The Poisson bracket (|J) is a particular case of the so- 
called Lie-Poisson bracket which differs from (||) by the 
presence of structure constants c a b c of some Lie algebra 
instead of e a b c characteristic of so(3). It is natural to 
think of the Nambu bracket as a particular case of 



{A,B,S} = c abc 



dA dB dS 

dx a dxb dx c 



(4) 
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where S = \g ab x a x is a second order Casimir invariant 
of an appropriate Lie algebra. Such brackets could be 
called the Lie-Nambu ones and are quite natural in the 
context of generalizations of a Lie-Poisson dynamics. It 
is surprising that this kind of generalization of Hamil- 
tonian dynamics has not been considered so far in the 
theory of classical dynamical systems (Ratiu, n.d.). One 
of the reasons seems to be the fact that for general Lie 
algebras the bracket does not satisfy the so called funda- 
mental identity (Chatterjee and Takhtajan, 1996). We 
shall return to this question in Sec. VII. 

The generalization of (g) which was extensively inves- 
tigated in the literature under the name of a "generalized 
Nambu dynamics" goes in another direction (Bayen and 
Flato, 1975; Takhtajan, 1994; Gautheron, 1996; Chat- 
terjee and Takhtajan, 1996). One treats the e not as 
structure constants but as a volume form. From this 
perspective it is natural to consider 



x H(b',b)p{c', 



fi, 



ibcH b p C 



(8) 



,A n } = 



dAi dA n 



(5) 



The parameter n is a dimension of the state space. It 
is not clear how to extend this type of description to 
infinite-dimensional spaces. 



B. Bialynicki-Birula— Morrison 1991: Liouville— von 
Neumann equation 

The observation that the Liouville von Neumann 
equation for a Wigner function can be written as a Lie- 
Nambu equation with nontrivial structure constants is 
due to Bialynicki-Birula and Morrison (1991). Below, 
instead of the Wigner function which is defined in terms 
of position-momentum coordinates, we shall stick to the 
more symmetric position-position representation. This 
will lead to a specific form of structure constants whose 
symmetry properties will be essential for further gener- 
alizations (Czachor, 1997a; Czachor and Kuna, 1997a). 

The density matrix in position representation is de- 
noted by p(a, a') =: p a , where we use a and a' instead of 
more typical x and x', and the lower composite index is 
introduced for brievity. The kinetic energy is represented 
by the kernel 



dyK(a,y)p(y,a') 



2m 



-p(a,a'), 



(6) 



and the Hamiltonian operator by 

H{a' 1 a)^K{a' 1 a) + V{a)5{a-a'). (7) 
It is easy to check that the following equation 

idtPa — J dbdb'dcdc' 

x(5(a- b')5{b - c')8(c - a') - 5{a ~ c')5(b - a')5{c - b') ) 



is equivalent to the Liouville- von Neumann one. The 
form (g) simultaneously illustrates the use of composite 
indices and the summation convention. Notice that the 
composite indices are in their lower or upper position, 
and the transition between the two is given by the metric 
tensor g ab working as follows 



pb = I dbdb'8(a — b')S(b — a')p(b, b') = p{a , a) = p a . 



!1 



So if p a = p(a,a') then p a = p(a',a). Although the lat- 
ter formula may seem somewhat artificial and was not 
used by Bialynicki-Birula and Morrison, it will prove ex- 
tremely useful when we arrive at various generalizations. 
The distributions Q a bc are structure constants of an 
infinte-dimensional Lie algebra, which can be checked by 
raising a with the help of g ab and verifying the standard 
properties. One should be aware of the fact that g ah is not 
the Cartan-Killing metric (which does not exist in this 
case). Writing g ab p a p = Tr(/5 2 ) —: C2 =■ 2S one rec- 
ognizes that g ab is a kernel form of the Hilbert-Schmidt 
metric. Let now H(p) = Tr Hp. Taking into account 
that H a = dH/8p a and p a = SS/8p a (S/5p a is a func- 
tional derivative) we can write the Liouville-von Neu- 
mann equation in the Bialynicki-Birula-Morrison form 
as 



tPk 



Pa wJ-£* s Ji = {p k ,H} 

OPb 0p c 



6p k SH SS 



{p k ,H,S}. 



(9) 

, , , u , (10) 

OPa OPb Opc 

Although the relationship of (g) and @ to (g) and (g) 
is obvious, it requires a few comments. First of all, the 
equations are Lie-Poisson and Lie-Nambu and not the 
"generalized Nambu" in the sense of the previous sub- 
section. The Lie algebra and the space of states are 
both infinite-dimensional. H is an average energy and 
C2 = 2S is a Casimir invariant called an "entropy" . 
(Actually, it is, up to constants, a 2-entropy of Daroczy 
(1970) or Tsallis (1988), and is closely related to Renyi's 
a-entropies (Renyi, 1960, 1961).) 



III. DIGRESSION ON PURE STATES 

There is one more formal prerequisite we need before 
we get further. The Liouville-von Neumann equation for 
mixed states has its roots in the pure-state Schrodinger 
equation. It turns out that the same is true of the struc- 
ture constants £l a bc- Let us switch now one level higher 
and instead of speaking about a nonrelativistic, spin-0 
Schrodinger equation consider general Hamilton equa- 
tions on a separable Hilbert space: 



SH 

Sip a > 



'h = |f. (id 
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IfV>a =ip(a),LJ aa = 6(a-a'),H = (ip\H\ip), and an ob- 
vious summation/integration convention is applied, (|ll| ) 
is equivalent to the Schrodinger equation (Chernoff and 
Marsden, 1974). ui aa is a symplectic form in the com- 
plex coordinates u ip = q + ip" . An explicit form of uj aa 
varies from representation to representation and is dif- 
ferent for, say, the Dirac equation, or a nonrelativistic 
particle with spin. It is important however that the form 
of the Hamilton equations (JTTJ) is always the same (al- 
though the dot at its LHS may have different meanings 
as well, cf. Czachor (1997a), Czachor and Kuna (1997a)) 
and that ui aa ip a (f> a > — (</>|'0)- (|ll| ) can be written in a 
form involving a Poisson tensor 

ilp a = Iaa'-^T-, -i$ a > = Icxa' TT"- ( 12 ) 
Oip a < Olj} a 



A pure-state density matrix is p a = p aa > = 4> a ^u' and 
u aa p a a' = u a p a = Tip. A pure-state Poisson bracket 
corresponding to ([l^) and its complex conjugated equa- 
tion is 



{A, B} = I a 



SA 6B 



8l/j a 8tp a 

SASB 
— Pall bc - , 

Opb 0p c 



{A «-> B) 



(13) 
(14) 



which holds for all functions A(p a ) — A(ip a ip a ') and 
B{pa) — B(ip a ?p a i). The structure constants are 



If we lower the indices in ( f20| ) we see that, somewhat 
counter-intuitively, g does not go directly into G (al- 
though g ab uj b = IJ) but 



bi...b n _ 



9a±...a n — Ga n a n _ 1 . , , ai . ( 22 ) 



So define a ^-operation which reverses the order of in- 
dices: ** = id, g* a = g a , g* ab = g ba = g ab , and 



9 ai...a n ~~ "fll»'<»n — 9a n a n -i...ai 



(23) 



and similarly with the upper indices. Essential for further 
calculations are the following properties. 

a) Cyclicity 



9&i ...a n b\. ,,b m — 9b\ . ..b m ai ...a n j 

* * 

y a\...a n b\...b m 9 b\...b rn a\...a rii - 



b) "Cut-and-glue" 



n 

ai-..a n x n a„ + i...a„ + „ 
y i)x 



c) "Annihilation'' 



J gai---akxa k+ i...a n + m _ ga 1 ...a k a k+ i...a n + r , 
W 9a 1 ...a k xa k+1 ...a, n - l - m 9a\ . . .a k a fc + i . . .a n +n 



d) "Drag-and-drop" 



(24) 
(25) 



(26) 
(27) 



(28) 
(29) 



Q a bc = 8f3 ,a Shifty 



where the deltas are defined by 



I a f3'UJ c 

7 , ,° 



= V Q 



(15) 
(16) 
(17) 



(18) 
(19) 



The metric tensor that raises and lowers the composite 
indices is g ab = w a/3 uj 13 " and g a b — I a p'Ipa'- The reader 
may check that we obtain the Bialynicki-Birula-Morrison 
formulas if wc replace cj's and Ps by the Dirac deltas. 



IV. HIGHER ORDER "METRIC" TENSORS 



gbi...bi ^ga 1 ...a k xa k + 1 ...a n+m _ g ai . . .a k bi . . .biak+i ■ ■ -a n+m ^gg-j 
9b\...bi 9a 1 ...a k xa k+1 ...a rl+m = 9a\ ...a k bi . . .bi a k + 1 . . .a„ +m ■ (31) 

It is practical to accept the rule stating that complex 
conjugation interchanges primed and unprimed indices. 
Assuming this we can define symmetric operators A as 
those whose kernels satisfy A a pr — Ap a i. We find also 
that 



3ai...a„ — 9 ai...a„ 

As a consequence 



(32) 



Qa-i.-.a^A^ 1 . . . A n n — Qai...a n A^ . . . A^ n (33) 

which is an abstract-index version of the well known rule 
Tr(i 1 ...l n ) = Tr(i n ...i 1 ) (34) 



In this Section we introduce several technical results 
which will turn our abstract composite index language 
into a practical tool. 



We have seen that Tr (p ) 



g PaPb- 



It is useful 



to introduce higher order tensors satisfying Tr(p n ) 
p an . Define 



ai...a„ 



valid for symmetric operators. In order to translate 
the abstract-index formulas into more standard operator 
ones one uses the following correspondence: 



a — .9aai...a„^i 1 ■ • • A n n 



(35) 



g ax...a n = yana^cja^osa, 
Gai...a n ^Qi a.' n ^-a.2Ct f , Iql^olL • 



(20) 
(21) 
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V. STRUCTURE CONSTANTS REVISITED 



VI. GENERALIZED LIE-NAMBU BRACKETS 



A Lie-Nambu 3-bracket written in the form (|4|) is based 
on a totally antisymmetric 3-index tensor. Obviously, the 
tensor has 3-indices for all Lie algebras and for this reason 
it is not immediately clear whether a generalization of (^) 
to a "generalized Nambu" n-bracket is possible. On the 
other hand, the structure constants occuring in dp have 
a rich structure and it turns out there exists a natural 
generalization of (|l0|). 

To begin with let us note that 



^abc 
n abc 



Qabc 
„abc 



Qacb : 
acb 



9 



2-9a[bc] 

2\g a ^ 



2 ] -9[abc], 

: 2\g^ ahc \ 



where [. . .] denotes an antisymmetrization. 
Consider 



fiai...a„ = ( n ~ l) ! S[ai...o„]- 



(36) 
(37) 



(38) 



Lemma 1. 



9[ ai ...a 2m ] - 
9[xai...a,2m] aa:[oi..43m]' 
^ X Qa 1 ...x...a n = 0. 



(39) 
(40) 
(41) 

Proof: (a) Eq. ©. 

Q[ai...a n ] 9[a-2 . . .a n ai\ ( 1) 9[ai...a rl ]i 

where the cyclicity and total antisymmetr y w ere used. 
The expression vanishes for even n. (b) Eq. (|0J). Assume 
n = 2m. 

9[xai ...a n ] 



2_t 

^ \9x[a 1 ...a„ 



•• • + ( 1) ' 9[a 1 ...a k \x\a k + 1 ...a„] + 

• • • + (-l) n fl[ Bl ...a B ]«) 

~ T ^ \9x[ai...a n ] + • • • + ( — 1) ' 9x[a k + i...a n ai...a k ] H~ • 

• • • H~ .9^c[ai . ,.a n ] J 

■ ■ • ~t~ 9x[ai*..a n ] J 



n + 1 



5x[(ii...fl Jl ] 

where we have used the cyclicity and the fact that (n — 
k + l)k is even for any k if n is even, (c) Eq. (f4l|). It is 
sufficient to note that the annihilation property together 
with @ and @ imply 



9[xai...a,2jn] 



Jx[a 1 ...a 2m ] — 9[a 1 ...a 2m ] 



□ 



We define the generalized Lie-Nambu bracket for n 
2m + 1 by 



{Ai, . . . , A ri \ — Hai...a n 

Let C k =g ai - ak p ai ...p ak . 
Theorem 1. 



SA X 5 An 

SPai ' ' ' SPa n 



(42) 



{Cfci, ■ ■ ■ 7 Cfc(„ + i) /2 ' • } — 0. 



(43) 



Proof: Let us begin with the following remark. Assume 
a tensor F... a b c ... has the "drag-and-drop" property 



F.. a 



be. 



C\C2 F abciC 



and consider 



F 



[o6c]...P 5 c lC2 P V 



g ( -fl ..!/(."■.. . "■ ' .../.'re/... ^ ...re//).. 



.acb... ~ F..coa... " F..oac...)pVc lC2 p C V C2 (44) 



F 



aci c^b. . 



F 



c\C2ba.. 



F.. baciC2 ...JpV 1 p C2 

- \(F... abciC2 ... - F. ClC2ba ...)pV 1 P C2 

which in general does not vanish. We can see therefore 
that for expressions such as (44) to vanish, it is sufficient 
to have in each term of the decompositition (Q) at least 
two transvected indices which are not seperated by a non- 
transvected one. With this observation in mind consider 



9x[a 1 ...a N c 1 ...c M ]9 a 1 ...a Kl 



■ <r 



P a i ... P a i 



(45) 



Expanding ( f45|) we will obtain a sum involving expres- 
sions g x ...ij... and g x ...ji... entering with opposite signs. 
If N and M will be chosen in a way guaranteeing that 
for any such term there exists at least a pair (q^. af) of 
indices that are not separated by some c r then (|5|) van- 
ishes on the basis of the preceding argument. We know 
that N + M is an even number. Therefore the great- 
est N that allows for a separation of any two cifc's is 
N = (N + M)/2 = M = (n - 1W2. To complete the 
proof it is sufficient to note that ( f45|) equals 



1 



'{Ck x ) • • • j Cfciv j '>•••) ' }> 



(46) 



where fe. 



Remarks: (1) Total antisymmetry of (E2L) guaran- 
tees that {Cfe,Cfc,...} = 0. Moreover (Ell) leads to 
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{Ci, . . .} = 0. Since C% — Ti p the dynamics generated 
by such brackets is trace preserving. (2) For any k 7 I 



{C k ,Q,...} 



pure states 



0. This follows from 



Q Pbx ■ • ' Pb n I pure states ' P @\ (pure states- 

(47) 

(3) Theorem 1 was proved for n = 3 in Czachor (1997a). 

Theorem 2. Let Sj — Sj(Ci, C2, ■ ■ ■) be a differentiable 
function of Ck, k = 1, 2, 3, . . . and z n a complex number. 
The dynamics given by 



Pa — z n {p a , Hi,...,H ( 



(n-l)/2j 



Si, 



1 S{n-X)/2) 



conserves Ch. Ch are Casimir invariants i.e. 



{C k ,A 1 



, ^(n-l)/2j 



Si, 



, S( n -l)/2} — 



(48) 
(49) 



for any functions A\ 



Remarks: (1) Theorem 2 is a straightforward conse- 
quence of Theorem 1. (2) The number z n will be as- 
sumed to satisfy z n = —z n (for n = Am + 3) or z n = z n 
(for n = Am + 1), m = 0, 1,2,.. .. The simplest choice is 
therefore either z n = —i, for n = Am + 3, or z n = 1, for 
n = Am + 1 (see the discussion below). 



VII. 3-BRACKET 

The simplest n = Am + 3 case is n — 3. The discussion 
given by Bialynicki-Birula and Morrison dealt with linear 
quantum mechanics. A possibility of using the 3-bracket 
dynamics as a departure point for nonlinear generaliza- 
tions of quantum mechanics was described in some detail 
in Czachor (1997a). One of the main motivations for 
studying the 3-bracket dynamics was a possibility of in- 
troducing nonlinearities only by generalizations of S and 
without modifications of H . Generalizations via nonlin- 
ear H are interesting and will be discussed in the next 
sections. An important drawback of such Hamiltonian 
generalizations is that we have to represent observables 
by nonlinear operators which leads to interpretational 
difficulties. To give an example, it is not clear which def- 
inition of a nonlinear eigenvalue is physically meaningful, 
or how to represent higher moments of experimentally 
measured random variables if nonlinear operators are in- 
volved (Czachor, 1996a). Let us therefore first consider 
what happens if H = H a p a = Tr Hp and S is an arbi- 
trary (differentiable) function of the Casimirs Ck- It is 
easy to see that the dynamics given by a 3-bracket is then 
linear if and only if S is linear in C2 (Czachor, 1997a). 



A. Lie-Nambu duality 

Is the 3-bracket a Poisson bracket? The answer to 
this question reveals an interesting duality which points 



into two different generalizations of linear quantum me- 
chanics. To understand the problem define {A, B}x '■= 
{A, B,X} and check whether the Jacobi identity is sat- 
isfied. Consider 

{{A,B} Xl C} x + {{C,A} X ,B} X + {{B,C} X , A} x 
5 A 5B S 2 X SC 5X 



Spd Sp e Sp a Spf Sp b Sp c 



(50) 



The terms involving second derivatives of A, B and C 
drop out just because of the total antisymmetry of struc- 
ture constants. The term involving the second drivative 
of X vanishes in several cases. For X = S = g ab p a p /2 
the second derivative gives g a f and ( |50| ) vanishes on the 
basis of the structure constants version of the Jacobi iden- 
tity. With this choice of X the bracket { • , • }s is a Lie- 
Poisson bracket and the dynamics given by 



Pa 



-i{p a , H} s = -i{p a , H, S} 



(51) 



is an ordinary Lie-Poisson dynamics. If H is nonlinear 
the dynamics corresponds to the nonlinear quantum me- 
chanics in the Bona- Jordan version (Bona, 1991; Jordan, 
1993). It can be shown that such brackets satisfy the 
Jacobi identity for all S = S(C 2 ) (Czachor, 1996b). If S 
is a function of higher order Casimirs, say, S = C3 , the 
Jacobi identity does not hold. However, rewriting (|51 



Pa = -i{p a , S}-h = -i{p a , S 7 -H} 



(52) 



we obtain a Poisson bracket for any S if X = —H is lin- 
ear. It follows that the requirement that observables are 
linear leads us, via the Lie-Nambu 3-bracket, to a dual 
Poisson structure given by ( |5^ ) . This Lie-Nambu duality 
(H, S) <-» (S 1 , —H) is typical of all Lie-Nambu theories 
and is analogous to the canonical (q,p) (p, — q) and 
electromagnetic (E, B) <-> (B, —E) dualities. The du- 
ality transformation is a particular case of the duality 
rotation 

{■, H,S} = { •, H cosa — S sin a, H sin a + iScosa}. (53) 

It is noteworthy that the since the A-bracket does not 
in general satisfy the Jacobi identity, the 3-bracket can- 
not satisfy the so-called fundamental identity discussed 
in Takhtajan (1994). The example of the Lie-Nambu 
duality, where the Jacobi identity simultaneously holds 
and does not hold (depending on the viewpoint), clearly 
shows that status of such identities is more technical than 
fundamental. Each kind of dynamics seems to have its 
own fundamental criteria of sensibility. In this work we 
insist on positivity of density matrices and lack of faster- 
than-light effects. 

Interesting in the context of the duality is the pure- 
state case where C„ = (ip\ip) n , S — S(\\ijj\\), and H = 
(tp\H\ip). The Hamilton equations (111) can be written as 
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itp a = H aa >—=^, -iipa' = H aa /——, (54) 
dtp a > dip a 

where / = uj aa ipa^a' = IIV'II 2 - If V = ^Pab is the 
electromagnetic spinor its squared norm \\ift\\ 2 is equal 
to the classical energy of the field. Therefore here H 
is not the energy although formally it is an expression 
analogous to the Dirac Hamiltonian function (Bialynicki- 
Birula, 1996). It is known that the Hamiltonian formal- 
ism based on the energy density ^(a;)! 2 = E{x) 2 + B(x) 2 
corresponds to the Poisson tensor which involves the 
differential operator J ■ iV = Vx (Bialynicki-Birula 
and Bialynicka-Birula, 1976). On the other hand tak- 
ing H = J ■ as the Hamiltonian function one 
gets the Poisson tensor which involves no differentiations. 
The nonlinear electrodynamics of the Born-Infeld type 
(Plebahski, 1970) may be regarded as a nonlinear gener- 
alization within the { • , S}-h scheme. 

Another manifestation of the duality can be seen in 
a simpler case of N harmonic oscillators. The classical 
energy E = J2(pI + €) = £ \Pk + «9fc| 2 = £ iV'fcl 2 can 
be regarded as a norm squared in the Hilbert space C N . 
It is an easy exercise to rewrite the equations of motion 
as an TV-dimensional Schrodinger equation with H being 
a diagonal matrix whose eigenvalues are the energies of 
the oscillators, but then H = (tp\H\ip) is not the energy. 

The above facts suggest an alternative interpretation 
of the Poisson structures that occur in quantum mechan- 
ics: A modification of H (say, by interactions) can be 
understood as a deformation of the Poissonian structure 
{ • , ■ }-h of the manifold of states, and not as a modifi- 
cation of a Hamiltonian function. Keeping H unchanged 
but modifying S one changes a flow on the Poisson mani- 
fold but the structure of the manifold itself is unchanged. 

B. Canonical transformations 



C. Formal solutions 

Consider first the dynamics with S = Ck+i/(k + 1). 
The 3-bracket equation 

Pa = -i{p a , H, S} = -in abc H b g cc ^- c *p Cl ...p Ck (55) 

has solutions which can be formally written as 

°° (— it) n 

= E —T {••■{{ A*(°)> H > 5 >> H,S}...,H, S}. (56) 

71=0 v ' 



Let p a (0) — p a (0). In the standard notation we have 
p(0) = p(0)* and 

'p = -i[H,p k \. (57) 

For p = p 2 ( |57j ) is the ordinary linear Liouville-von Neu- 
mann equation. Asuming H* = H we find p(0)* = p(0). 
In the same way we can prove that d n p/dt n \ t= o = 
(d n p/dt n \ t= o)* . It follows that the formal solution sat- 
isfies p(t) — p{t)* if p(0) = p(0)*. The same argument 
applies to more general S = S(C\, C*2, ■ ■ ■)■ ( |57| ) is inter- 
esting in itself even in finite-dimensional cases where the 
above argument can be made more rigorous. To show 
that spectrum of self-adjoint Hilbert-Schmidt solutions 
of the 3-bracket equations of motion is conserved by the 
3-bracket dynamics one uses the following 

Lemma 2. Consider a sequence of probabilities {pk}kLo 
and an arbitrary real sequence {cik}%° =0 satisfying for any 
natural n 

oo oo 

= (58) 

k=Q k=0 

Then the two sequences are identical up to permutation. 

Remarks: (1) Let the two sequences represent spectra 
of a Hermitian Hilbert-Schmidt solution of the n-bracket 
equation at t — and t ^ respectively. Since t \— > p(t) is 
continuous, the spectrum of p(t) is also continuous and, 
hence, conserved. The condition (|5^) is implied by con- 
servation of C n . (2) Lemma 2 was proved in Czachor and 
Marciniak (1997). (3) Such solutions can be interpreted 
as nonlinear ly evolving density matrices. The question 
of their complete positivity will be discussed below. 



Nonlinear quantum mechanics based on { • , • }<j uses 
nonlinear H and S = S^C^). The canonical transfor- 
mations must therefore be those that do not change 
C-2 = Tr(/5 2 ) (or (ip\4>) for pure states). Such transfor- 
mations can be nonlinear and were discussed by Wein- 
berg (1989) and Jordan (1993). The version based on 
{ • , • }—h leads to canonical transformations that keep 
H = TrHp linear (or (ip\H\ip) bilinear for pure states). 
The two classes of transformations are not equivalent. It 
is natural to require that only those observables which 
commute with H have to be represented by linear func- 
tionals. Various Poissonian structures that appear in 
this context may be used also for a combined quantum- 
classical description as shown by Jones (1993, 1994) for 
the Weinberg type theory. 



VIII. 5-BRACKET 

n = 5 is the simplest (nontrivial) 4m + 1 case. The 
equation of motion is (z§ = 1) 

p a = {p a ,Hi,H 2 ,Si,S 2 }. (59) 

Assume that are linear in p. The simplest choice 
of the other two generators is S\ — Ca/2, S-2 = C3/3. 
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Nontrivial 5-bracket equations of motion (^) are always 
nonlinear, as opposed to the 3-bracket ones which can be 
linear, and always vanish on pure states. The RHS of 
( |59l ) when written in the standard notation involves an 
antisymmetrizcd product of Hi, H 2 , p and p 2 . After a 
simplification one finds 

'p={\p,Hi]H 2 -\p,H 2 ]Hi)p 2 
+p 2 (H 2 [Hi 1 p]-Hi[H 2l p\) 
+p(H 2 p 2 Hi - Hip 2 H 2 ) 



+ (Hip 2 H 2 -H 2 p 2 H 



IP- 



(60) 



The RHS of (|6^) is Hermitian if p, Hi and H 2 are Her- 
mitian. This explains the choice of real z§. For p = p 2 
( |60| ) vanishes. Assuming that p(0) = p(0)* we find that 
all higher derivatives are also Hermitian. The formal so- 
lution 

p«(*) = X)-r 

n=0 

x{...{{ Pa (0),Hi,H 2 ,Si,S 2 },Hi,H 2 ,Si,S 2 } 

n 

...,Hi,H 2 ,S u S 2 } (61) 



satisfies p(t) = p(t)* if p(0) = p(0)*. Using the same 
argument as for n — 3 we conclude that the spectrum 
of self-adjoint and Hilbert-Schmidt solutions of (|59|) is 
conserved. 



IX. iV-PARTICLE EXTENSIONS OF 1-PARTICLE 
ALMOST-LIE-POISSON DYNAMICS 

An extension of dynamics from 1 to AT particles is a 
delicate problem. Careful analysis shows that the Lie- 
Nambu duality holding for the 3-brackets leads to gen- 
eralizations which behave differently from the viewpoint 
of A-particle extensions. The Poisson dynamics based 
on { ■ , • }s for S = C 2 /2 is the most regular one. An 
inclusion of nonlinear Hamiltonian functions H does not 
lead to difficulties with independent evolutions of sepa- 
rated systems. This fact was proved by Polchinski (1991) 
and Jordan (1993) and recently generalized by myself to 
those nonlinear theories which do not possess Hamilto- 
nian functions but only Hamiltonian operators (Czachor, 
1997b). In spite of this, the view that any nonlinear gen- 
eralization of a Schrodinger dynamics leads to problems 
with causality is quite popular. Nonlinear quantum me- 
chanics based on { • , -}-h leads to a new kind of nonlocal 
phenomenon. This effect, typical of mixed states, is anal- 
ogous to the threshold phenomena discussed by Goldin 
and Svetlichny (1994) for pure states. 



A. iV-particle "metric" tensors 

Let g ai ...a n be a 1-particle metric tensor. The N- 
particle tensor is defined by 



-.N 



9a}. ..oi ■ ■ ■9a?...a»- 



(62) 



The indices at the LHS of (|62J) are the iV-particle ones: 
di = a\ . . . ai , etc. The TV-particle n-bracket is defined 
by 

(A A \N C>N SA l SA n 

{Ai,...,A n } =n ai _ anT - w ... T - w , (63) 



Spg Sp* 



where p£ = p a i... a * 



^...a n = (n-l)!<,.. a „] 



(64) 



For identical particles (bosons and fermions) p a i,.. a N = 
P(a 1 ...a N ) with (. . .) denoting symmctrization. A distinc- 
tion between fermions and bosons can be seen at the 
"spinor index" level: 

Pa 1 . ..a" = P[a 1 \a' 1 \...\a N ]a' N = Pa 1 [a' 1 1 . . . \a N \a' N ] (fermions), 
Pa 1 . ..a" = P(a 1 \a' 1 \...\a N )a' N = Pa 1 (a' 1 1 . . . \a N \a' N ) (boSOns) . 

For the 3-bracket the annihilation property implies the 
important identity 



. . . <jj - - <jj ■■ ' - . . . OJ 
— 9aiC! ■ ■ ■ 9a k ^ 1 c k _ 1 Qa k b k c k ga k + 1 c k+1 ■ ■ -9a N c N , (65) 

where £l ak b k c k are the 1-particle structure constants of 
the fc-th particle. 



B. Extension of Hamiltonians 

The 3-bracket dynamics is Lie- Poisson if S = C 2 /2 and 
a Hamiltonian function H exists. There exist nonlinear 
Liouville-von Neumann equations that possess Hamilto- 
nian operators of the form H(p) = Hi(p) + H 2 (p) where 
TiH 2 {p)p = 0. Such equations do not possess a Hamil- 
tonian function Tr H(p)p but often are of physical inter- 
est (e.g. starting with non-lincarizable Doebner-Goldin 
Schrodinger equations (Doebner and Goldin, 1996) one 
arrives at this class of mixed state equations (Czachor, 
1997b)). Equations that can be written as 



p a = -i£labcH b (p)p c 



(66) 



although no H satisfying H b = |^ exists, will be called 
almost-Lie-Poisson. Assume we have N (not necessarily 
identical) particles that do not interact with one another 
(but can interact with something else and do not have to 
be free). Each of them satisfies a 1-particle equation (69) 
with some H. We define the iV-particle extension of (6q) 
by 
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p N a 



(67) 



where 
H\ P N ) 



= H^ 1 (p (1) V b2 . . . u b « + . . . + ^ . . . Lo b »-iH b N » (p w ). (68) 
The reduced density matrix p;u\ is defined by 



P(k)a k = U ai 



...uj aN p ai ... ak ... aN . (69) 



There are two motivations for (pq) . First, if the fc-th 
particle is described by a Hamiltonian function Hk{p) = 
Hk{P(k)) then (|68|) is just a consequence of the chain rule 
for functional derivatives. The second motivation is (65). 
Indeed, applying (|65|) to (p7|) we obtain 



if>a 1 ...a N - Qa 1 b 1 c 1 H bl (p( 1 ))p Cl a2 ...a N 

+ ■■■ + n aNbNCN H b N »(p {N) )p ai ... aN _ 1 CN . (70) 



Transvecting both sides of (f70J) with 



,a k _ 1 a k+1 



and using uj a Vl a i, c = we get 

P(k)a = -itta k b k c k H bk (P(k))P 



Cfc 

(fe)- 



(71) 



Both sides (|7l]) depend only on objects which are intrinsic 
to the fc-th subsystem. It follows that the reduced den- 
sity matrix of this subsystem "does not see" the other 
noninteracting systems. The observers in the other sub- 
systems have no possibility of influencing the dynamics 
of the k-th one by any kind of modification of the Hamil- 
tonians in the other separated systems. In particular, 
they cannot influence any observable quantity in the k-th 
subsystem by different choices of measurements in their 
"own" subsystems. This explicitly contradicts the popu- 
lar claim that any nonlinear dynamics must imply faster 
than light influences between separated systems. 

Denoting the dynamics of the iV-particle system by 
(f)^, the one corresponding to the k-th subsystem by (j) k , 
and by Tr jv_fc the partial trace which reduces the dy- 
namics from the composite system to the fc-th one, we 
get an important separability condition 



Tr N _ k o^ = ^oTr N _ k 



(72) 



characteristic of the Lic-Poisson dynamics. The dynam- 
ics satisfying (jT^) and (|7l]) can be termed completely sep- 
arable. 



C. Examples of completely separable extensions 



The method of extension given by (|68|) applies to 
any equation whose 1-particle Hamiltonian operator can 
be written as a function of the particle's density ma- 
trix. This applies also to pure state (Schrodinger) equa- 
tions. To see how this works for nonlinear Schrodinger 



equations consider the examples. A Hamiltonian oper- 
ator consists of two parts: The linear part Hl(x) — 
-^kinetic + V(x) and a nonlinear Hnl = H(ip,tp;x). To 
apply the above method we have to be able to write 
H(ip, 4>] x) as H(p; x). 

a) "Nonlinear Schrodinger equation" . 

H(ip,-ip;x) = \ip(x)\ 2 -> H(p;x) = p(x,x) 

b) Bialynicki-Birula-Mycielski equation (Bialynicki- 
Birula-Mycielski, 1976). 

H(ip, "tjj; x) — In (|^(x)| 2 ) — > H(p; x) = lnp(x, x) 

Obviously in the same way one can treat any equation 
with nonlinearities given by some function H(\ip(x)\). 

c) Haag-Bannier equation (Haag and Bannier, 1978) 

tp(x)W x ip(x) - ip(x)V ' x i>(x) 



H{i>,i>;x)=A{x) 



2i\^(x)\< 



= A W 2ip~fax) (?3) 

d) Doebner-Goldin equations (Doebner and Goldin, 1996; 
Nattermann, 1997). There are five nonlinear terms de- 
noted by Rk- 



Ri(p;x) 



1 ip(x)A x ip(x) - ip{x)A x ip{x) _ 
2i \^(x)\ 2 ' 

1 J dzS(x — z)A x [p(x, z) — p(z, xj\ 
2i p(x,x) 



D , ,, 7 \ &x\^(x)\ 2 A x p(x,x) 
i? 2 (V>, ip; x) = L) / ^ |2 > R 2 {p; x) = 



R 3 (i/),ip;x) 



\m\ 2 

i [^(x)v x ^(x) - i>(x)v x i>(x)] ; 



(2if 



\Mx)\ 



1 (J dz8(x- z)V x [p(x,z) -p(z,x)]) 



Ri^ylp-yX) 



(2z) 2 p{x,x) 2 
1 [^{x)V x ^{x)-ip{x)V x 4>{x)} -V x \tP(x)\ 2 _ 
2i IV'O)! 4 
1 J dzS(x - z)V x [p(x, z) - p(z, x)] ■ V x p(x, x) 
2i p{x,x) 2 



, - \\lM(x)\ 2 \^ , x \V x p(x,x)Y 
R 5 (^;x)= 1 , ,\ \ J ^R 5 (p;x)- - 



\ip(x)\ 4 ^' ' Pix.xf 

e) Twarock equation on S 1 (Twarock, 1997) 



H(ip,ip;x) 
H(p;x) 



ip(x)"ip(x)' - i)(x)"ip(x)' 
ip(x)ip(x)' — ip(x)ip(x)' 

[J dyS(x - y)d 2 p(x, y)} [J dz5(x - z)d 2 p{z, x)] - c.c. 
p{x, x) J dy5(x - y)d x p(x, y) - c.c. 
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f) (n, n)-homogeneous nonlinearities . Denote by D a 
differential operator involving arbitrary mixed partial 
derivatives up to order k. Consider a real function 
H (tp) — F[Di/j(x)j , which is (n, n)-homogeneous i.e. sat- 
isfies H(Xtp) = X n \ n H(tp). We first write 



F(Dip(x)) 



F(ip(x)Dip(x)) 
\ip{x)\ 2n 



and then apply the tricks used for the Haag-Bannier, 
Doebncr-Goldin and Twarock terms. Obviously any rea- 
sonable function of such (n, n)-homogeneous expressions 
with different n's will do as well. 

Let us now concentrate on the simplest case with 
Hk{p;x) = Hk(p(x,x)) and just two particles. The 2- 
particle extension of the nonlinear part of the Hamilto- 
nian is 

+ H 2 (p (2) (x 2 ,x 2 )) 

= Hi(Jdyp(xi,y ) x 1 ,y) + H 2 (jdyp(y,x 2 ,y,x 2 ). (74) 

If the 2-particle state is pure p(x\, x 2 , x\, x' 2 ) = 
^(xi,x 2 )^*(x' 1 ,x' 2 ) the RHS becomes 

Hi(/dy|*(a!i,tf)| 2 ) + H 2 {fdy\*(y,x 2 )\ 2 ) (75) 

and reduces to 



H^ix^l 2 ) + H 2 (\4>(x 2 )\ 2 ) 



(76) 



on product states &(xi,x 2 ) — ip(xi)(j)(x 2 ). 

It seems that an example that cannot be treated in 
this way is the Kostin equation (Kostin, 1972) involving 
nonlinearity In (ip{x) /ip{x)) . 

A reader may have noticed that the above reasoning 
involves two "heretic" elements. First of all the 2-particle 
extension of dynamics for nonfactorizable (entangled) 
states leads to integro-differential equations. Such equa- 
tions are typically rejected in nonlinear quantum mechan- 
ics literature as nonlocal. The construction presented 
above shows that it is in fact just the opposite. The 
requirement of locality (complete separability) leads us 
to appropriate integral terms and precisely because of 
these terms the subsystems can be completely isolated 
from one another. Second, all nonlinearities of the form 
F(\ip(x)\ 2 ) are acceptable. This is in an apparent con- 
tradiction with the well known result of Bialynicki-Birula 
and Mycielski who used the separability criterion to de- 
rive the logarithmic nonlinearity. However, they assumed 
that the 2-particle extension has to be F(|^(xi, x 2 )\ 2 ) 
and, with the condition 

F(|V>(zi)#r 2 )| 2 ) = F(|V(xi)| 2 ) + F(\^(x 2 )\ 2 ), (77) 

they found that only F ~ In is acceptable. One of the 
obvious drawback of such extensions is that they do not 
tell us what to do if the systems are noninteracting but 
correlated and *5>(xi,x 2 ) does not factorize. In such a 
case local probability densities are obtained by integrat- 
ing out the coordinates of the remaining particles and 
it is quite logical that such expressions occur in the N- 
particle "correctly extended" Hamiltonians we discussed. 



D. Problem of complete positivity 

A subsystem described by p(M can be embedded into 
a composite one described by p N in a way guaranteeing 
a consistency of (^) and ([n]). The dynamics of p/Q is 
independent of N. In addition, since both and p N 
satisfy the 3-bracket Lie-Poisson equation, the extension 
procedure preserves positivity of dynamics at both sub- 
system and composite system levels. A dynamics that 
has these properties is typically associated with the no- 
tion of a completely positive map, provided the maps are 
linear. 

In mathematical literature the notion of complete pos- 
itivity is generalized to nonlinear maps in a way that can 
be translated to our context as follows (Ando and Choi, 
1986; Arveson, 1987; Majewski, 1990; Alicki and Majew- 
ski, 1990). One takes a positive map 



4>\{a) = a(t), 4>\ : A -> A 



(78) 



where A is a unital C*-algebra. In our case a — p^) 
and 4>\ (a) — p^){t). Assume for simplicity that the di- 
mension of the fc-th system is finite. In the next step 
one considers a density matrix p N (0) of a bigger system 
consisting of the original one plus a system which has a 
finite number m of degrees of freedom. Writing 



P N (0)= E P N (0)rT'ss'\r}(r'\®\s)(s'\ 

r,r' ,s,s' 

= ^o ss - ® \s)(s'\ 



we can represent p N (0) by the matrix 
( an • ■ ■ ai„ 

V ^ml ■ ■ • (Lmi 



(79) 



(80) 



whose entries are elements of A. The (nonlinear) map 4>\ 
is said to be completely positive if the matrix 



/ 0*i ( a n) ... 4>\ (aim) 

\ 0l(«ml) ■ ■ • 4>\{ a rnm) 



is positive for any m. This is equivalent to the positivity 
of 



^(p>))=^i(a ss 0®| S >< S '| 



(82) 



However, for nonlinear cf>\ the explicit form of ( p2| ) for 
t > is different for differenet choices of bases {|s)}, 
which is unphysical unless there exists a superselection 
rule distinguishing a particular basis. In the generic case 
no such distinguished basis exists. Therefore a basis inde- 
pendent extension from 1 to more particles cannot have 
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the forms (^lj) and (]82)- And, indeed, the dynamics fol- 
lowing from the Lie-Poisson extension discussed above 
does not coincide with fl82]). This was shown by an ex- 
plicit calculation in Czachor and Kuna (1997b) but could 
be inferred also from the basis-independence of the N- 
particle extension. It must be stressed that the dynamics 
( |S2| ) is the one that was used by Gisin (1989) in his discus- 
sion of unphysical influences between separated systems. 

X. NONLOCAL PROPERTIES OF iV-PARTICLE 
EXTENSIONS FOR THE DUAL POISSON 
STRUCTURE {•, -}-h 

The regularity of the TV-particle extensions typical of 
an almost-Lie-Poisson dynamics is lost when one consid- 
ers the dual Poisson structure { • , S}-h with S being a 
higher order Casimir invariant. To explicitly see the kind 
of difficulties one may encounter consider the 2-particle 
equation 

ip ai a 2 — {p ai a 2 , Cs/3}-H- (83) 

General properties of the 3-bracket dynamics imply that 
C n (p) are conserved for any natural n, where p is the 
2-particle density matrix. Also Ci{piW) is a constant of 
motion. However 

iC 2 {p {1) ) =2Tr 1 ([Tr 2 (p 2 ),Tr 2 (/3)] J ff 1 ), (84) 

where the indices 1 and 2 correspond to the two sub- 
systems and we have assumed the standard 2-particle 
extension of the (linear) Hamiltonian. Although we do 
not have much control over the behavior of the eigen- 
values pj of the reduced density matrix pm, we can 
infer that ^2jPj is constant whereas YljPj i s > m S en ~ 
eral, time dependent. Let us note that average energies 
of the two subsystems are separately conserved. This 
follows from the general property of the 3-bracket: For 
H(p) = H 1 {p {1) )+H 2 (p (2) ) 

{H 1 (p (1) ),H(p),S(p)} = (85) 

for any S (Czachor, 1997a). Therefore the probabilities 
Pj can be made time dependent without making the two 
subsystems interact with each other and without chang- 
ing energies of the subsystems just by modifying the over- 
all entropy of the composite system. So the change of 
entropy, say, by C 2 — » C 2 + eC*3 at the global level leads 
to the modification of the local subsystems. Such a mod- 
ification will not occur if 

[Tr 2 (/5 2 ),Tr 2 (p)] =0 (86) 

which holds for a pure-state p, or p — pn\ ® /3( 2 ). Still, 
strong correlations can also help since reduced density 
matrices occuring in a singlet state are proportional to 
unit matrices and the commutator vanishes. Systems de- 
scribed by entropies other than C 2 /2 possess some kind 



of overall identity which is lost when it is physically 
meaningful to discuss their subsystems separately. This 
effect deserves a name. The fact that the subsystems 
"feel" that the total entropy (information) undegoes a 
change from C 2 to C 2 + eCs although apparently "noth- 
ing happened" (no energy has been transfered between 
the neighboring subsystems) resembles the influence the 
Big Brother from G. Orwell's "1984" exerted on the in- 
habitants of Oceania "just by watching them" . 

It is not immediatly clear that the "Big Brother effect" 
is entirely unphysical. Its interpretation is obscured by 
our lack of understanding of the physical role played by 
the entropies C n in the multiple-bracket scheme. It may 
be relevant to mention that C 2 is characteristic of Renyi 
2-entropy which is the only a-entropy that characterizes 
a system whose gain of information is for all proba- 
bility distributions. The analysis of this problem was 
given by A. Renyi (1960). Although this is the first pa- 
per where the notion of a-entropies was introduced, it 
does not seem to be known to the majority of experts 
in quantum mechanical information theory. The work 
typically quoted in the literature is Renyi (1961) where 
Renyi already departed from the natural definition of the 
information gain in favor of a "decrease of uncertainty" . 
This latter modification was motivated by the problem 
with the vanishing gain for a — 2. 

A class of physical systems whose identity as a whole 
is associated with the way their entropy (or information) 
behaves are living organisms. Similarly, statistical prop- 
erties of societies have dynamical properties strongly de- 
pending on information, and their dynamics cannot be 
regarded as a simple sum of individual activities. The 
fact that a possibility of gaining information can be for- 
mally related, via S, to nonlinearity of evolution resem- 
bles a similar phenomenon mentioned by Wigner (1967) 
in the context of the measurement problem ("paradox 
of a friend" ) . Whether such phenomena are in any way 
related to the 3-bracket dynamics is at the moment a 
matter of pure speculation. 

XL QUANTIZATION OF CLASSICAL NAMBU 
DYNAMICS? 

The (2n -I- l)-bracket can be regarded as a nonlin- 
ear quantization of a classical (n + l)-bracket with n 
classical Hamiltonian functions H%, . . . , H n . Indeed, the 
Liouville-von Neumann equation is characterized by one 
Hamiltonian operator Hi, obtained by a quantization of 
a classical Hamiltonian function H i . The requirement of 
linearity of evolution combined with the 3-bracket dy- 
namics leads to the choice of S = C 2 /2. Having n 
Hamiltonian functions Hk we can obtain n Hamiltonian 
operators Hk after some quantization procedure (say, 
p — » — iKVx, etc.). Representing the operators by kernels 
Hi*, k = 1, . . . , n, we can consider the (2n + l)-bracket 
equation 
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Pa — Z2n+lQ 



■ Hi 



Tja, 



5S n 



5Sx 

Spbi Spb n 



(87) 



For n > 1 the equation is always nonlinear and its RHS 
vanishes on pure states. A self-adjoint Hilbert-Schmidt 
solution of the "quantized Nambu dynamics" may be in- 
terpreted as a density matrix because spectrum of the 
solution is conserved. Had we started with the linear 
equation, which could be obtained by taking the (rt + 2)- 
bracket 



Pa 



z n+2 



H n P 



(88) 



we would have obtained a dynamics which would not, in 
general, conserve Tr (p m ) for m > 1 and there would be 
no guarantee that positivity of p is conserved. 

This kind of nonlinear quantization differs from the 
procedure discussed in Takhtajan (1994) which was based 
on an n-bracket obtained by an antisymmetrization of a 
product of n oparators, or the Zariski product quantiza- 
tion proposed in Ditto et al. (1997). Also all operator 
expressions involving an odd number of operator kernels, 
if described within our approach, must be excluded be- 
cause the "metric" tensor used for the generalized struc- 
ture constants would have to have an even number of 
indices, but such structure constants vanish (3-bracket 
involves antisymmetrization of two operators, 5-bracket 
antisymmetrizes four operators, etc.). The quantization 
proposed originally by Nambu (1973) (cf. Garcia Sucre 
and Kalnay, 1975) is therefore also incquivalent to our 
formulation. 



XII. WHAT NEXT? 

The formalism presented in this work is at a very pre- 
liminary stage of development. The main problem is 
how to solve the nonlinear density matrix equtions and 
how to extend the approach to a fully relativistic the- 
ory. Both questions are highly nontrivial. The equations 
of the form ip — [H , p n ] bear some formal similarity to 
the Nahm equations studied in the SU(2) monopole the- 
ory (Hitchin, 1983). Some recently developed techniques 
of solving matrix equations by a noncommutative ver- 
sion of a Darboux transform (Leble and Zaitsev, 1997; 
Leble, 1997) may prove useful in this context. A natural 
candidate for a relativistic multiple-bracket formalism is 
the off-shell proper-time formulation. Some preliminary 
indefinite-metric results can be found in Czachor and 
Kuna (1997a), and a work on positive-metric generaliza- 
tion of the Bargmann-Wigner off-shell equations (Cza- 
chor, 1997c) is in progress. The characteristic function 
off-shell approach developed recently by Naudts (1998) 
seems especially suited for this kind of generalization. 
A separate problem is the behavior of eigenvalues of re- 
duced density matrices in the { • , -}-h scheme. These 
eigenvalues are in general time dependent and therefore 
the question of their interpretation is still unclear. 
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